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The lifetime of the Bs meson is measured using the semileptonic decay Bg 
Djt^uX. The data sample consists of 19.3 pb~^ of pp collisions at i/s = 1.8 TeV 

collected by the CDF detector at the Fermilab Tevatron collider during 1992-1993. 
There are 76 it 8 i~^D~ signal events where the Dg is identified via the decay D~ 
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^7r~, (f) — K^K" . Using these events, the Bg meson hfetime is determined to be 
Ts = 1.42 to 23 (stat) lb 0.11 (syst) ps. A measurement of the Bs hfetime in a low 
statistics sample of exclusive Bs — > J/'ip4' decays is also presented in this paper. 
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The lifetime differences between the bottom hadrons can probe the -B-decay mechanisms 
which are beyond the simple quark spectator model. In the case of charm mesons, such 
differences have been observed to be quite large {t{D^)/t{D^) ~ 2.5). Among bottom 
hadrons, the lifetime differences are expected to be smaller due to the heavier bottom quark 
mass. Phenomenological models predict a 5-10% difference between the and meson 
lifetimes and very similar Bd and Bs lifetimes [||]. This is consistent with the previous 
measurements of Bu^d meson lifetime f^, as well as recent Bg lifetime measurements from 
LEP 1^ . It has also been suggested by recent theory calculations |Q that the lifetime between 
the two CP eigenstates produced by mixing of the Bg and Bs may differ by as much as 15%. 
Such an effect may manifest itself as a difference in lifetimes between the Bg semileptonic 
decay, which is almost an equal mixture of the two CP states, and the decay Bg J/'4"Pi 
which is expected to be dominated by the CP even state. In this letter, we first present 
the measurement of Bg lifetime using the semileptonic decay p Bg Dji'^uX, where the 
is identified via ^tt^, </> K^K^. We then describe briefly a result using the 

exclusive decay Bg J/ip'Pi where J/iIj ^ f^^fJ' , <P —>■ K^K". The data sample for this 
paper consists of 19.3 pb'^ of pp collisions at y/s=1.8 TeV collected by the CDF detector 
during the 1992-1993 run. 

The CDF detector is described in detail elsewhere [|] . We describe here only the detector 
features most relevant to this analysis. Two devices inside the 1.4 T solenoid are used for 
the tracking of charged particles: the silicon vertex detector (SVX) and the central tracking 
chamber (CTC). The SVX consists of four layers of silicon microstrip detectors located at 
radii between 3.0 and 7.9 cm from the interaction point and provides spatial measurements 
in the r-(p plane [0 with a resolution of 13 /im, giving a track impact parameter resolution of 
about {13 + AO /pt) /im [^], where px is the transverse momentum of the track in GeV/c. The 
transverse profile of the beam is circular and has an RMS of ~ 35 /im, while the longitudinal 
beam size is ~ 30 cm. The CTC is a cylindrical drift chamber containing 84 layers grouped 
into 8 alternating superlayers of axial and stereo wires. It covers the pseudorapidity interval 
\ri\ < 1.1, where rj = — ln[tan(^^/2)]. The Pt resolution of the CTC combined with the SVX 
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is S(pt)/pt = ((0.0066)2 + (0.0009pt)^)^/^- Qutside the solenoid are electromagnetic (CEM) 
and hadronic (CHA) calorimeters (|?7| < 1.1) that employ a projective tower geometry. A 
layer of proportional wire chambers (CES) is located near shower maximum in the CEM and 
provides a measurement of electromagnetic shower profiles in both the if and z directions. 
Two different muon subsystems in the central region are used, the central muon chambers 
(CMU) and the central upgrade muon chambers (CMP), with total coverage of 80% for 
|?7| < 0.6. The CMP chambers are located behind 8 interaction lengths of material. 

Events containing semileptonic Bg decays were collected using inclusive electron and 
muon triggers. The Et threshold for the principal single electron trigger was 9 GeV, where 
Et = Esin{6) and E is the electromagnetic energy measured in the calorimeter. The single 
muon trigger required a pt > 7.5 GeV/c track in the CTC with matched track segments in 
both the CMU and CMP systems. 

Offline identification of an electron involved measurements from both the calorimeters 
and the CTC. Photon conversion electrons were removed by searching for an oppositely 
charged track which had a small opening angle with a primary electron candidate. 

A muon candidate was required to be detected by both the CMU and CMP chambers 
to reduce background due to hadrons that do not interact in the calorimeter. Good po- 



sition matching [|Ty] was required between track segments in the muon chambers and an 
extrapolated CTC track. 

The 07r~ reconstruction started with a search for candidates. We first defined a 

search cone around the lepton candidate with a radius AR = y/Arp~+~A(p^ of 0.8. Any two 
oppositely charged tracks with pt > 1 GeV/c within that cone were assigned kaon masses 
and combined to form a candidate. No kaon identification was used in the selection. 
Each candidate was required to have Pt{4>) > 2.0 GeV/c and a mass within ±8 MeV/c^ 
of the world average mass [|TT| . The candidate was then combined with another track of 
Pt > 0.8 GeV/c inside the cone which had the opposite charge of the lepton (the 'right-sign' 
combination). This third track was assigned the pion mass. To ensure a good decay vertex 
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measurement, track quality cuts were imposed on the lepton and at least two of the three 
track candidates forming the Dg candidate. The ^ K~ , and tt" tracks were then refit 
with a common vertex constraint. The confidence level of that fit was required to be greater 
than 1%. Since the has spin 1 and both the D~ and tt~ are spin 0, the helicity angle 

which is the angle between the and directions in the rest frame exhibits a 
distribution dN/dlcos'^) ~ cos^ \E'. A cut |cos\I'| > 0.4 was therefore applied to suppress 
the combinatorial background, which we found to be a flat in cos \1/ distribution. The mass 
of the iDs system was required to be between 3.0 and 5.7 GeV/c^ in order to be consistent 
with coming from a Bg decay. We also applied an isolation cut E^°/pT{4>'n'^) < 1.2 on the 
candidate, where E}^° is a sum of transverse energy within a cone of radius 0.4 in r]-ip 
space around the lepton candidate, excluding the lepton energy. This cut eliminated many of 
the fake combinations from high track multiplicity jets. Furthermore, we required that 
the apparent Dj decay vertex (VdJ be positively displaced from the primary vertex along 
the direction of the i^Dj momentum. Figure |I|a shows the (f)TT~ invariant mass distribution 
for the 'right-sign' lepton-Dg combinations. A Dg signal with mean of 1.967 GeV/c^ and 
width of 5.4 MeV/c^ is observed. Evidence of the Cabibbo suppressed D~ (pn^ decay is 
also present. No enhancement is seen in the corresponding distribution for the 'wrong-sign' 
combinations (Figure Pd). We select a signal sample using a D~ mass window of 1.953 to 
1.981 GeV/c^. A total of 139 events are found with a background fraction f^g = 0.45 ±0.01. 
The number of (-^D~ events above background in the sample is 76 ± 8. 

There are two possible sources of non-strange B meson decays which can lead to right- 
sign C.^D~ combinations. The first one is a four body decay Bu,d D^'Ki'^u, where K 
denotes any type of strange meson. Because of the low probability of producing ss pairs and 
the limited phase space, this process is suppressed and has not been observed experimentally. 
The recent ARGUS limit (90% CL) is BR(5„,d D-Kf+u) < 1.2% p[. Also, a theoretical 



analysis based on the 'resonance model' yields BIi{Bu,d ~^ D^K.i^u) < 0.025 x BK{Bd — »• 
i'^uX) |T^. Using the latter result and our estimated efficiency, we expect less than 2.6% of 
our i^D^ combinations from this source. The second process is B^^d ~^ D^DX, D —>■ i'^uX, 
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where D is any charmed meson. This decay produces softer and less isolated leptons than 
that from Bg semileptonic decay and therefore the acceptance for this source relative to 
the signal is quite small (~ 2.6%). Using the BR(i? DgX) |1^,[15| and the semileptonic 



branching ratios of and [|TT|, we estimate the fraction of this type of background 



is less than 3%. In addition, we also considered the background from cc production where 
a DjD pair is produced. Monte Carlo predicts the background fraction from this type of 
source to be < 7%. In summary, the contribution of all above physics backgrounds is quite 
small compared to the combinatorial background. We will consider them as a source of 
systematic uncertainty for the Bg lifetime measurement. 

The secondary vertex where the Bg decays to a lepton and a D~ (referred to as VbJ is 
obtained by intersecting the trajectory of the lepton track with the flight path of the D~ 
candidate. The transverse decay length L is defined as the displacement in the transverse 
plane of Vb^ from the primary vertex projected onto the direction of the pxiiDs). This 
is our best estimator of the Bg direction. The effect of the unknown Bg relativistic boost 
can be partially removed event-by-event with the factor pT{iDg)/M{Bg) (where M{Bg) = 
5.37 GeV/c^ [0) and leads to a corrected decay length 

^ PTiiDg) ' 

which is referred to as the 'proper decay length'. A residual correction between pxiiDg) and 
Pt{Bs) is done statistically by convoluting a Monte Carlo distribution of the px correction 
factor K= pt{(!-D g) / pt{B g) with an exponential decay distribution in the lifetime fit. The K 
distribution has an average value of 0.86 and an RMS of 0.11 and is approximately constant 
as a function of pxi^Dg). To model the proper decay length distribution of the background 
events contained in the signal sample, we define a background sample which consists of 
the right-sign events from the D~ sidebands (1.885-1.945 and 1.990-2.050 GeV/c^) and the 
wrong-sign events from the interval 1.885-2.050 GeV/c^. 

The proper decay length distribution (Figure 0) is fit using an unbinned maximum log- 
likelihood method. Both the Bg lifetime and the background shape are determined in a 



simultaneous fit using the signal and background samples. Thus the likelihood function L 
is a combination of two parts: 

Ns Nb 

>c = n[(i - M^s^g + hgHg] ■ n (2) 

j i 
where Ns and Nb are the number of events in the signal and background samples. The 
signal probability function J^sig consists of a normalized decay exponential function (defined 
for only positive decay lengths and symbolized by £^4.) convoluted with the K distribution 
and a Gaussian resolution function Q: 

^sigic-T. s) = £+{-Kx, ct) ® K''"* ® Q{C - X, sa'), (3) 

where ^* is the measured proper decay length with uncertainty a* (typically lOO/im) and 
X is the true proper decay length. The scale factor s accounts for the underestimation of 
the decay length error. The background is parameterized by a Gaussian centered at zero, 
symmetrical positive and negative exponential tails, and a positive decay exponential to 
characterize the heavy flavor decay in the background sample. 

The best fit values of ct and s are found to be 426 ^ gg /zm and 1.4 ± 0.1 respectively. 
Figure ^a shows the proper decay length distribution of the signal sample with the result 
of the fit superimposed. The same distribution of the background samples is shown in 
Figure 03. As a consistency check, we also fit the Dg lifetime from the proper decay length 
measured from the tertiary vertex Vd^ to the secondary vertex Vb^- The result is ct{Ds) = 
135 1 

30 which is consistent with the world average value [|ll|. 
Table I lists all sources of systematic uncertainty considered in this analysis. Major 
contributions come from the source of the background shape, the non-i^^ production, and 
the resolution function. To model the contribution to the signal from the combinatorial 
backgrounds, we combined the events from three different sideband regions. There may be 
some bias in choosing the correct mixture. We find a ± 4% variation in the lifetime when 
using each sideband region individually. The dominant source of systematics from non-i?^ 
production was found to be B^^ D^DX decays. This mode was studied using Monte 

10 



Carlo simulations and the contribution to the systematic uncertainty in the lifetime was 
found to ± 4%. The effect of the decay length resolution was studied by varying the scale 
factor and using an alternative resolution function consisting of two Gaussians, giving a 3% 
systematic uncertainty. 

Quoting the statistical and systematic uncertainties separately, we measure the Bg life- 
time using semileptonic decays to be 

TB, = 1.42 t'oil (Stat) toll (syst) ps. 

This result is consistent with the previous world average of 1.34 IIqIj ps |[TT[| . 

For the exclusive mode measurement, we use the decay chain Bg J/'ip(py J/i' ^ 1^^ i 
(f) K^K^ . The data sample and reconstruction techniques used for this decay channel 
are similar to those described in detail elsewhere |T^. Briefly, the invariant mass of two 



oppositely charged muon candidates is calculated after the tracks are constrained to originate 
from a common vertex. J /if) candidates are selected by requiring the difference between 
the dimuon mass and the world average J/ip mass |TT|] to be < 3(T, where a is the mass 
uncertainty calculated for each dimuon candidate. The meson selections are the same as 
reference but with Pt{<P) > 3.0 ( rather than 2 GeV/c) to further reduce the background. 
To reconstruct Bg meson candidates, the 4 daughter tracks are constrained to originate from 
a common vertex and the dimuon mass is simultaneously constrained to the world average 
J/ip mass. We require the probability for this combined fit to be > 2%. In addition, at 
least one of the /i candidates and at least one of the other tracks must be well measured in 
the SVX. 

Using the measured pt{Bs), the proper decay length is calculated and the lifetime of 
the Bg meson is determined by performing a simultaneous unbinned log-likelihood fit to the 
entire mass and proper decay length spectra. The mass distribution is fit to a Gaussian and 
a flat background. We model the proper decay length distribution of the background with 
a Gaussian centered at zero and positive and negative exponential functions. The signal is 
described by an exponential decay function convoluted with a Gaussian. This fit determines 
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the mass and lifetime of the Bg, the signal fraction, and background shape simultaneously. 
The proper decay length distribution is shown in Figure ^, where we have displayed events 
within ±21 MeV/c^ of the Bg mass peak. The fit returns 7.9 i'^ g signal events and a Bg 
hfetime of ctb, = 520 taio /^m. 

We estimate a total systematic error of ±20/im, the dominant contributions arising from 
the uncertainty in the parametrization of the background shape and our understanding of 
the resolution function. The Bg lifetime using fully reconstructed Bg J/'4"P decays is 
measured to be: 

tb, = 1.74 tosl (stat.) ± 0.07 (syst.) ps. 

In conclusion, the Bg lifetime has been measured in both the semileptonic and exclusive 
decay channels. At present, the two measurements are consistent with each other within 
their quoted uncertainties and are consistent with the results of the Bu and Bd lifetimes 
previously measured by CDF. 

We anticipate a more significant result in both modes after the ongoing collider run. 

We thank the Fermilab staff and the technical staffs of the participating institutions for 
their vital contributions. This work was supported by the U.S. Department of Energy and 
National Science Foundation; the Italian Istituto Nazionale di Fisica Nucleare; the Ministry 
of Education, Science and Culture of Japan; the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council of the Republic of China; the A. P. Sloan 
Foundation; and the Alexander von Humboldt-Stiftung. 
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TABLES 



TABLE I. Semileptonic mode systematic uncertainties. 



Systematic Source 


Uncertainty 


Background shape 


4% 


Non-^s source 


4% 


Resolution function 


3% 


Boost correction 


2% 


Decay length cut 


2% 


Fitting method 


1% 


Misalignment 


2% 


Total 


7% 
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FIGURES 




1.8 1.9 2 2.1 2.2 



M(^7Y ) [GeV/c'] 

FIG. 1. 

The mass distribution of 07r^ for (a) 'right-sign' combination {(f)TT~£^); (b) 'wrong-sign' combina- 
tion {(j)7r~ i~).The shaded regions are used for the background sample. 
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FIG. 2. 

(a) Proper decay length distribution for the t^Dj signal sample with a curve (solid) from unbinned 
log-likelihood fitting of signal and background. The dashed line represents the contribution from 
the combinatorial background. The dotted one represents the signal contribution, (b) The proper 
decay length distribution for the background sample with a curve representing the background 
lifetime contribution. 
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FIG. 3. 

Bg lifetime measurement using the J/ipcf) exclusive mode. Inset figure is the mass distribution. The 
solid curves show the fit results. 
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